Wannier-Stark states of a quantum particle in 2D lattices 
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A simple method of calculating the Wannier-Stark resonances in 2D lattices is suggested. Using 
this method we calculate the complex Wannier-Stark spectrum for a non-separable 2D potential 
realized in optical lattices and analyze its general structure. The dependence of the lifetime of 
Wannier-Stark states on the direction of the static field (relative to the crystallographic axis of the 
lattice) is briefly discussed. 



I. INTRODUCTION 

The quantum states of a particle in a periodic po- 
tential plus homogeneous field (known nowadays as the 
Wannier-Stark states, WS-states in what follows) are one 
of the long-standing problems of single-particle quantum 
mechanics. The beginning of the study of this problem 
dates back to the paper by Bloch of 1929, followed by 
contributions of Zcner, Landau, Wannier, Zak and many 
others [GJ. In the late eighties the problem got a new 
impact by the invention of semiconductor superlattices. 
The unambiguous observation of the WS-spectrum in a 
semiconductor superlattice Q ended a long theoretical 
debate about the nature of WS-states, and now it is com- 
monly accepted that they are the resonance states of the 
system. Besides, WS-states were recently studied in a 
system of cold atoms in an optical lattice || and some 
other (quasi) one-dimensional systems. 

Although WS-states are resonances, i.e. metastable 
states, in the theoretical analysis of related problems 
they were usually approximated by stationary states 
(one-band, tight-binding, and similar approximations). 
Beyond the one-band approximation, WS-states in the 
semiconductor and optical lattices were studied in recent 
papers W and || by using the scattering matrix approach 
of Ref. || (see also Ref . for details) . This approach ac- 
tually solves the one-dimensional Wannier-Stark problem 
and supplies exhaustive information about ID WS-states. 
In the present letter we extend the method of Ref. I], |?|] 
to the case of two-dimensional lattices. For the first time 
we find the complex spectrum of 2D WS-states and ana- 
lyze its general structure. 

To be concrete, we choose the following system: 



H = p 2 /2 + V(r) + F • r 



0,2/) 



V(r) = cos x + cos y — e cos x cos y , 



(1) 
(2) 





where < e < 1 joj. Two limiting cases e = and 
e = 1 correspond to an 'egg crate' potential, for which 
the system is separable, and a 'quantum well' potential, 
where the coupling between two degrees of freedom is 
maximal (see Fig. |]). Let us also note that the choice 
e = 1 corresponds to a 2D optical potential created by 
two standing laser waves crossing at right angle. Thus 
the results presented below can be directly applied to the 
system of cold atoms in a 2D optical lattice. 



II. 2D WANNIER-BLOCH SPECTRUM 

We briefly recall the key points of the ID theory. The 
spectrum of the Bloch particle in the presence of a static 
field consists of several sets of equidistant levels 



£ a ,i =E a + 2nFl - iT a /2 



(3) 



known as Wannier-Stark ladders of resonances. In 
Eq. (||), 2-7T stands for the lattice period, F is the am- 
plitude of the static force, I = 0, ±1, ... is the site index 
and the index a — 0, 1, . . . labels different ladders. The 
lifetime of WS-states ^ a ,i(x) is defined by the resonance 
width r a as r Q = h/T a . Typically, the lifetime r Q rapidly 
decreases with increasing index a. Because of this only 
the first few WS-ladders are of physical importance. 

Along with the WS-states ^ at i(x), one can also intro- 
duce Wannier-Bloch states (WB-states) by 



ip a ,k(x) = 22^g t i(x) exp(i27rfcQ 



(4) 



FIG. 1: Potential energy m for e = (a) and e = 1 (b). 



As follows from the definition (|j), the continuous evo- 
lution of WB-states obeys the equation ij) a ^(x,t) — 
exp(-i£ a t/h)ip aik -Ft/h(x), where £ a = E a - iT a /2. 
Thus, WB-states can be alternatively denned as the 
eigenfunction of the evolution operator over the Bloch 
period Tg = K/F j|. (Note that the eigenvalues of the 
evolution operator form degenerate bands £ a (k) = £ a ). 
Additionally, to ensure that i/> a ,k(%) are resonance states 
of the system, the eigenvalue equation for the evolution 
operator should be accomplished by the specific non- 
hermitian boundary condition. It was proven in Ref. R] 
that the required boundary conditions are imposed by 
the truncation of the evolution operator matrix in the 
momentum representation. 
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We proceed with the two-dimensional case. As men- 
tioned above, WB-states in a ID lattice can be defined 
as the non-hermitian eigenstates of the evolution opera- 
tor over one Bloch period. In the 2D problem there are 
two different Bloch periods associated with the two com- 
ponents of the static field. Therefore the notion of the 
WB-states can be introduced only in the case of com- 
mensurate periods, i.e., in the case of 'rational' direction 
of the field (q,r are coprime integers): 



F x = 



qF 



( r 2 +q 2y/2 ' 



rF 



( r 2 +g 2)l/2 



(5) 



Provided condition (||) is satisfied, we define 2D WB- 
states as the non-hermitian eigcnfunctions of the sys- 
tem evolution operator over the common Bloch period 
Tb = i 7 ' 2 + q 2 ) 1 ^ 2 ^/ F . Using the Kramers-Henneberger 
transformation, which is just the gauge which transforms 
the static term into the vector potential, the evolution 
operator can be presented in the form 



U(T B ) 



H(t) = 



-iix e -ir» gjg 



(PX - F X t) 2 (fiy 



dtH(t) 



(6) 



Fytf 



+ V(x,y), (7) 



which reveals its translational invariance (the hat over 
the exponent sign denotes time ordering). Alternatively, 
we can rotate the coordinates so that the direction of the 
field coincides with the x'-axis: 



qx + ry 
(r 2 +q 2 ) 1/2 ' 



y' = 



qy 



( r 2 +(? 2)l/2 



(8) 



Transformation (ph introduces a new lattice period a = 
2-7r(r 2 + q 2 ) 1 / 2 and reduces the size of the original Bril- 
louin zone s = r 2 + q 2 times. Associated with the new 
lattice period is a new Bloch time T a = [r 2 +q 2 )~ 1 / 2 h/F, 
which is s times shorter than the original Bloch time Tg. 
Using p' x = —ihd/dx' and p' y = —ihd/dy', the time evo- 
lution operator over the new Bloch time T a in the rotated 
coordinates has the form 



U'(T a ) = e- i27nr '/ Q exp 



dtH'(t) 



(9) 



H'(t) 



(p' x Ft) 2 _ + v^ + v{x ,^^ 



2 2 
Then, presenting the wave function as 



Mr') = e lkV V c n , (r'ln') , (r'ln') = - e i2 ™'- r > 



we get the matrix equation 

V[/ ,(k,) c , - e- i£Ta ' h c , 



(11) 

(12) 



where U^J, denotes the fc'-dependent matrix elements 
of the operator : 



U'¥] = (n' 



,-ik'-: 



U'(T a )e 



ik'- 



m 



(13) 



Similar to the ID case, the truncation of the infinite uni- 
tary matrix (fL3j), 



K\,\<\<N- 



- oo 



n' 



\m' y \ < M 



■ CO 



(14) 



which is presumed in the numerical calculations, auto- 
matically imposes the non-hermitian boundary condition 
along the x'-direction. (Truncation of the matrix over the 
index n' y , m' y does not change the hermitian boundary 
condition along the y'-direction.) Then the eigenvalues 
£ obtained by numerical diagonalization of the truncated 
matrix correspond to the quantum resonances. 

In the transformed coordinates, the unit cell with area 
a 2 = (2ir) 2 s contains s different sublattices, and each of 
them supports its own WB-states. The sublattices are 
related by primitive translations of the unrotated lattice, 
and correspondingly the energies of their WB-states dif- 
fer by multiples of aF/s. Furthermore, as function of 

the quasimomentum, the energies £ — £p (k' x , k') (here 
j3 = 0, 1, ... is the 'Bloch band' index and i = 1, . . . ,s is 
the sublattice index) do not depend on k^. This follows 
from the fact that a change of k' x in Eq. ([13|) can be com- 
pensated by shifting the time origin in Eq. (^). For the 
y'-degree of freedom the Bloch theorem can be applied, 
and therefore £^ ' {k' xl k' y ) is a periodic function of k' y with 
generally nonzero amplitude A£g. Thus, assuming a ra- 
tional direction of the field, in each fundamental energy 
interval aF, the static field induces s = r 2 + q 2 iden- 
tical sub-bands, separated by the energy interval aF/s. 
Simultaneously, the size of the Brillouin zone is reduced 
by a factor s. This result resembles the one obtained 
for a ID lattice affected by a time-periodic perturbation 
JlOt or that for a 2D lattice in a magnetic field [^l). In 
these cases - provided the condition of comensurability 
between the Bloch period and the period of the driv- 
ing force or the condition of 'rationality' for the mag- 
netic flux through a unit cell, respectively, is fulfilled - 
the (quasi)energy spectrum of the system has a similar 
structure. 

We conclude this section with a remark concerning the 
numerical procedure. Although the reduced Brillouin 
zone approach described above is the most consistent, 
we found it more convenient to diagonalize the evolu- 
tion operator without preliminary rotation of the coordi- 
nate. In other words, in order to find the WB-spectrum, 
we solve the eigenvalue equation (|l2| ) with the truncated 
matrix constructed on the basis of the operator (||). As 
a result of the diagonalization, one obtains eigenvalues 
£p{kx,k y ) with quasimomentum k = (k x ,k y ) defined 
in the original Brillouin zone. Because the WB-bands 
are uniform along the direction of the field, £p{k x ,k y ) 
is a periodic function of both k x and k y with periods 
1/r and 1/q respectively. The energies obtained in this 
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FIG. 2: Position of the ground WB-band repeated by the 
subband energy interval 2-kF(t 1 + g 2 )~ 1//2 as a function of 
the field direction 9 = arctan(r/g) (parameters Ti — 2, F — 
0.08^/2, t = 0, integers q,r < 21). 
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FIG. 3: Real (left) and imaginary (right) parts of the dis- 
persion relation £p(k x , k y ) for the ground WB-states and dif- 
ferent values of the potential parameter e = 0, 0.1,0.5, and 
1 (from top to bottom). The system parameters are Ti — 2, 
F x = F y = 0.08, and k y = 0. 



way can then be used to construct the complete WB- 
spectrum £ g\k' x , k' y ), i = l,...,s. In the next section 
we present results of a numerical calculation of the dis- 
persion relation £p(k x ,k y ) for the periodic potential (0) 
and moderate values of the static field F = (F x ,Fy), 
|F| = F = const. 

III. NUMERICAL RESULTS 

It is instructive to begin with the separable case e = 0. 
In this case, 2D WB-states are given by the product of 
ID states and 2D WB-energies are just the sum of ID 
energies. In what follows we restrict ourselves to analyz- 
ing only the ground band. First we consider the real part 
of the spectrum Eq — Re(£rj)- 

It was shown in the previous section that for rational 
directions of the field the ground WB-subbands repeat 
with energy splitting aF j s. As an example, Fig. || shows 
the relative positions of these subbands as a function of 
the angle 8 = &rct&n(r/q) for h = 2 and F = 0.08^2. We 
recall that in the considered case of a separable potential 
the bands have zero width for any 8 ^ 0, w/2. 

The main difference between separable and non- 
separable potentials is that the subbands Eq\1l) have 
a finite width in the latter case. This is illustrated by 
Fig. ||(a) which shows the dispersion relation Eo(k x , k y = 
0) for the potential (||) with (from top to bottom) e = 0, 

0. 1, 0.5, and 1. The direction of the field is 8 = ir/4, 

1. e. r = q = 1. The amplitude of the static field and 
the value of the scaled Planck constant are the same as 
in Fig. ||. It is seen in Fig. |(a) that the WB-bands 
gain a finite width as e is increased. We also calculated 
the dispersion relation Eo(k x ,k y = 0) for different an- 
gles 8 = arctan(r/<7), with r, q < 6. It was found that the 



band widths AEq = AEq (r, q) are typically much smaller 
than the mean energy separation between the subbands. 
Thus, for practical purpose, one can neglect the band 
width for the real part of the spectrum. (An exception is 
the case 8 = 0, 7r/2 where the width of the WB-bands ap- 
proximately coincides with the width of the Bloch band 
in the absence of the static field.) Neglecting the width 
of the bands they were found to form a structure similar 
to that shown in Fig. ||. 

We proceed with the analysis of the decay rate of the 
WB-states, which is determined by the imaginary part 
of the complex energy, To = — 2Im(£o)- In the case of 
a separable potential the dependence Tq = To (F, 8) is 
obviously given by the equation 

T (F,8) =T' (Fcos8)+T' (Fsm8) , (15) 

where T' (F') stands for the width of ID WS-resonances. 
For the parameters used (K = 2 and F — 0.08V2) the 
dependence ([l5]) is shown in Fig. [| by a solid line. The 
maximum around 8 = Tt/2 originates from a peak-like 
behavior of T' (F') and is explained by the phenomenon 
of ID resonant tunneling [Q. 

For a non-separable potential and rational direction 
of the field the decay rate depends on the quasimomen- 
tum. For the particular case 8 — 7r/4 this dependence is 
depicted in Fig. ||(b). We would like to note the compli- 
cated behavior of To(k). The oscillating character of the 
decay rate is an open problem for the present day. Be- 
cause the decay rate depends on the quasimomentum it 
might be convenient to introduce the notion of Po, where 
the average is taken over the reduced Brillouin zone. The 
dots in Fig. ^ show the values of f o for some rational di- 
rection of the field and two different values of e. It is 
seen that for a small e = 0.1 the ratio AFo/Fo is small 
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FIG. 4: Decay rate of the ground WB-states as a function of 
the field direction 6 in the case of separable potential (e = 
0, solid curve). The dashed and dashed-dotted lines are an 
interpolation to arbitrary d of the mean decay rate calculated 
for some rational directions of the field (dots) for e = 0.1 and 
e = 1, respectively. The maximum and minimum values of 
the decay rate for these angles are indicated by the 'error' 
bars. 

and the obtained dependence Tq = Tq(t, q) essentially 
reproduces that of the separable case. However, this is 
not valid for e = 1, where the decay rate varies wildly. 



Thus, in the case of strong coupling between two degrees 
of freedom the description of WS-state by a mean decay 
rate is insufficient. 



IV. CONCLUSION 

We studied Wannier resonances in a 2D system, mainly 
discussing the complex energy spectrum of the Wannier- 
Bloch states. However, because the latter are related to 
the Wannier-Stark states by a Fourier transformation, 
the obtained results can be easily reformulated in terms 
of the Wannier-Stark resonances. Then the following is 
valid, (i) Neglecting the asymptotic tail, WS-states are 
localized functions along the direction of the field. (This 
follows from the degeneracy of WB-bands along the field 
direction.) (ii) For any rational direction of the field 
[see Eq. (g)] WS-states are Bloch waves in the transverse 
direction, (iii) For a non-separable potential the corre- 
sponding energy bands have a finite width, (iv) For the 
real part of the spectrum, the band widths are small and 
can be well neglected for r, q > 1. 

We also found a nontrivial dependence of the resonance 
width (inverse lifetime of WS-states) on the direction of 
the field. Because the value of the resonance width de- 
fines the decay of the probability, a complicated behavior 
of the survival probability is expected when the direction 
of the field is varied. The detailed study of the probabil- 
ity dynamics is reserved for future publication. 
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